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ABSTRACT 

We present new results on the kinematics, thermal and ionization state, and spatial distribution of 
metal-enriched gas in the circumgalactic medium (CGM) of massive galaxies at redshift ~ 3, using 
the "Eris" suite of cosmological hydrodynamic "zoom-in" simulations. The reference run adopts a 
blastwave scheme for supernova feedback that produces large-scale galactic outflows, a star formation 
recipe based on a high gas density threshold, metal-dependent radiative cooling, and a model for the 
diffusion of metals and thermal energy. The effect of the local UV radiation field is added in post- 
processing. The CGM (defined as all gas at R > 0.2R v { r = 10 kpc, where R v { r is the virial radius) 
contains multiple phases having a wide range of physical conditions, with more than half of its heavy 
elements locked in a warm- hot component at T > 10 5 K. Synthetic spectra, generated by drawing 
sightlines through the CGM, produce interstellar absorption line strengths of Lya, C n, C iv, Si n, 
and Si iv as a function of galactocentric impact parame ter (scaled to the v irial radius) that are in 
broad agreement with those observed at high-redshift by iSteidel et al.l (|2Q1Q[ ) . The covering factor of 
absorbing material declines less rapidly with impact parameter for Lya and C iv compared to C n, 
Si iv, and Si n, with Lya remaining strong (WL ya > 300 mA) to > 5R v i Y = 250 kpc. Only about one 
third of all the gas within R v [ r is outflowing. The fraction of sightlines within one virial radius that 
intercep t optica l ly thi ck, TVhi > 10 17,2 cm -2 material is 27%, in agreement with recent observations 
bv iRudie et al.l (j2012[ ). Such optically thick absorption is shown to trace inflowing "cold" streams 
that penetrate deep inside the virial radius. The streams, enriched to metallicities above 0.01 solar 
by previous episodes of star formation in the main host and in nearby dwarfs, give origin to strong 
(iVcn > 10 13 cm -2 ) C ii absorption with a covering factor of 22% within R v [ r and 10% within 2R Ylv . 
Galactic outflows do not cause any substantial suppression of the cold accretion mode. The central 
galaxy is surrounded by a large O vi halo, with a typical column density A^ovi ^ 10 14 cm -2 and a near 
unity covering factor maintained all the wa y out to 150 kpc. This ma tches the trends recently observed 
in star- forming galaxies at low redshift bv iTumlinson et al. I (j2011f ). Our zoom-in simulations of this 
single system appear then to reproduce quantitatively the complex baryonic processes that determine 
the exchange of matter, energy, and metals between galaxies and their surroundings. 
Subject headings: galaxies: evolution - galaxies: high-redshift - intergalactic medium - method: nu- 
merical 



1. INTRODUCTION 

Studies of the ionization, chemical, thermodynamic, 
and kinematic state of gaseous material in the circum- 
galactic medium (CGM) hold clues to understanding the 
exchange of mass, metals, and energy between galax- 
ies and their surroundings. It is the poorly understood 
flows of gas into galaxies, and from them back into their 
environments, that determine the response of baryons 
to dark matter potential wells, regulate star formation, 
and enrich the intergalactic medium (IGM). Observa- 
tions clearly show that galactic-scale metal-rich outflows 
with velocities of several hundred kms -1 are ubiquitous 
in massive star- for ming: galaxies at high redshift (e.g., 
iPettini et all 1200 H iShaplev et all 120031 : IVeilleux et all 
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120051 : IWeiner et all I2009D . Far-UV spectra of close 
angular pairs of z ~ 2-3 "Lyman Break Galaxies" 
(LBGs) have recently provided a detailed map of such 
metal- enriched gas as a f unction of galactocentric dis- 
tance (jSteidel et al.l [2010). and a benchmark for sim- 
ulations of the galaxy-IGM ecosystem. The observed 
line strengths of several ionic species convey information 
about the physical conditions, covering fraction, and ve- 
locity spre ad of the absorbing multi-phase CGM on 100 
kpc scales (Steidel et al. 2010). These, together with the 
increasing amount of data that has been accumulated 
recently o n the galaxy environs at high and low red- 
shifts (e.g. Prochaska & Hennawi 2009; Crighton et al. 
201 It [Pr ochaska et al. I 120111 : ITumlinson et al. I 1201 H 



Rudie et al. 2012), must be reproduced by hydrodynam- 



ical simulations of galaxy formation that aim at follow- 
ing the transport of heavy elements from their produc- 
tion sites into the IGM (e.g. iQppenheimer fc D ave 2008; 
Wiersma et all 12001 IShen et al.l 120101 : ICen fc Chisaril 



20lTh . 



Recent theoretical work has found that the major- 
ity of baryons accreted by galaxies below a certain 
critical mass flow in "cold" rather than being shocked 
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to the virial tem p eratu re (e.g. |B irnboim fc Dekell 

20031: lKereJ_el_aLJ |2QQl [Dekel & Birnboim 2006; 
Qcvirk. Pichon. fc Tevssierl 120081 : IDekel et all [2009; 



Keres et al. 1 120091 ). This cold accretion mode occurs 
mainly along dense streams that deliver fresh fuel fo r 
star formation to the disk (e.g. iBrooks et al. I [2009). 
Even when a shock is present, cold gas accretion can 
occur along fila ments that penetr ate dee p inside the 
hot halo (e.g. IDekel fc Birnboiml 120061 : iDekel et all 
120091 lAgertz. Tevssier. fc Moord 120091 ). Signatures of 
such streams may include Lv a blobs dGoerdt e t al. I 
12010( 1. Lyman Limit Systems (iFu magalli lFaD 120111 ). 
or Mg ii absorbers (jKacprzak et al. I l2012f ). Previ- 
ous numerical studies on cold flows and the cycle of 
baryons in and out of galaxies have either focused 
on a statistical description of how galaxies get their 
gas and sacrificed n umerical resolution on individual 
galaxy scales (e.g. |Keres et al. I 120051: lOcyirk et al. I 
2008; IDekel et al-Tl2QQ9b Ivan de Voort et al. 1 12011( 1. or 
used "zoom-in" simulations that do not generate strong 
galactic winds and hence cannot address the complex 
interaction between infl owing and outflowing gas in the 
multiphase CGM (e.g . iFaucher-G iguere & Keres 20111 
i Fumagalli et al. 11201 lUKimm et al. 112011b Stewart et al 
I2OIII: iGoerdt et al. Il20l3v ! 

In IShen et aTl (j2Q12h , we presented a high-resolution 
analysis of the metal-enriched CGM of a massive galaxy 
at z = 3 using results from the "Eris" suite of new cosmo- 
logical hydrodynmic "zoom-in" simulations, in which a 
close a nalog of a Milky Way system arises at the present 
epoch (jGuedes et al.ll20llh . Such simulations adopt a 
blastwave scheme for supernova feedback that - in com- 
bination with a star formation recipe based on a high 
gas density threshold - gener ates galactic outflow s with- 
out explicit wind particles. She n et al.l (|2012f ) identi- 
fied three sources of heavy elements within 200 physical 
kpc from the galaxy's center: the main host, its satel- 
lite progenitors (accreted by the main host before red- 
shift 3), and nearby dwarfs. The CGM was shown to 
be a blend of outflowing metal-rich and infalling metal- 
poor gas, with approximately half of all gas-phase met- 
als locked in a hot- warm component. In this Paper we 
use a new simulation of the "Eris" suite to create syn- 
thetic far-UV absorption spectra through the CGM and 
compare them to high-redshift observations, and to make 
predictions about the properties and detectability of cold 
accretion flows. To anticipate the results of our study, 
the simulated transmission spectra are found to be in 
good agreement with the interstellar metal absorption 
line strengths and H 1 covering factors as a function of 
impact parameter recent ly me asured in a z =2-3 galax y 
sample bv iSteidel etaLl (j2010l ) and iRudie et al.l ()2Q12h . 
We show that cold accretion streams within 1-2 R v [ r are 
traced by optically thick, metal-enriched absorbers with 
covering factor of 10-20%. There is no substantial sup- 
pression of the cold accretion mode caused by galactic 
outflows. The central galaxy is surrounded by a large 
halo of collisionally- ionized O vi plasma that matches the 
properties of the O vi bubbles recentl y detected around 
l ow-re dshift starf-forming galaxies by iTumlinson et al. I 
(|20lTh . 

2. SIMULATION 



i: 



The Eris suite of zoom-in simulations of a Milky 
Way galaxy analog is being performed in a Wilkin- 
son Microwave Anisotropy Probe 3-year c osmology with 



the p arallel TreeSPH code Gasoline (Wa dslev et al 
l2Q04f) . Details of main run were given in llGuedes et al 
(j2011f ). The high-resolution region, 1 Mpc on a side 
at z = 0, contains 13 million dark matter particles 
and an equal number of gas particles, for a final dark 
and gas particle mass of tubm — 9.8 x 10 4 M and 
^sph = 2 x 10 4 M , respectively. The spline soft- 
ening length is 120 pc. The twin run analyzed here 
("Eris2") presents severa l improvements over the previ- 
ously published versions (jGuedes et al.ll20lTI : IShen et aTl 
120 121 ): 1) the initial stellar m ass function (IM F) follows 
the modern determination bv iKroupal (j2001f ). This in- 
creases the number of Type II supernovae per unit stel- 
lar mass by about a factor of 2, and the IMF- averaged 
metal yield by a factor of 3 compared to iKroupa et al.l 
(1993); 2) the metallicity-dependent radiative cooling 
at all temperatures in the range 100- 10 9 K is deter- 
mined using pre-computed t abulated rates from the pho- 
toi onization code Cl oudy (jFerland et al.lll998l ). follow- 
ing |Shen~eF_2] (|201Q[ ). Cloudy tables assume that met- 
als are in ionization equilibrium. The ionization, cool- 
ing, and heating rates for primordial species (H, H + , 
He, He + , He ++ ) are calculated time-dependently from 
the rate equations. A uniform, cosmic TJV background 
modifies the ionization and excitation state of the gas, 
photoionizing away abundant metal ions and reducing 
the cooling efficiency. I t is implemented using the new 
lHaardt fc Madaul (j2Q12h redshift-dependent spectra, in- 
cluding emission from quasars and star- forming galaxies. 
A non-uniform local stellar radiation field is added in 
post-processing, as detailed below. As in previous runs, 
the gas is assumed to be optically thin to ionizing ra- 
diation at all wavelengths; 3) star formation occurs sto- 
castically when cold (T < 10 4 K), virialized gas reaches 
a threshold density of usf = 20 atoms cm -3 (4 times 
higher than used in the original Eris run), and proceeds 
at a rate dp*/dt = 0.1p gas /£dyn, where p* and p gas are the 
stellar and gas densities, and td yn is the local dynamical 
time; 4) star particles inject energy, mass, and metals 
back into the ISM through Type la and Type II SNe and 
stellar winds. We track the formation of Oxygen and Iron 
separately, and convert Oxygen to alpha-elements and 
Iron to i ron-peak elements a ssuming solar abundances 
patters (jAsplund et al.l [20 09). In the absence of some 
implementation of diffusion, metals are locked into spe- 
cific parti cles and their distribu tion is artificially inhomo- 
geneous (Wi ersma et al.l 12009). Eris2 includes a scheme 
for turbulent mixing in shearing flows that redistributes 
heavy elements and thermal energy b etween win d mate - 
rial and the ambient gas. Following IShen et al.l ((2010), 
the mixing of any scalar quantity A is approximated by 
a diffusion term 



dA 
~dt 



V(DVA), 



(1) 



D 



where th e diffusion co efficient D follows the model pro- 
posed by Smagorinskyj (|1963f ) for the atmospheric bound- 
ary layer, 



D 



C \Sij\ h\ 



(2) 
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Fig. 1. — The growth of Eris2's metal enriched CGM. The figure shows the projected gas metallicity (i.e. the surface mass density of 
metals divided by the total mass surface density) at different redshifts in a cube of 600 (proper) kpc on a side. The galaxy stellar disk is 
seen nearly edge-on in this projection. At z = 2.8, the galaxy center is indicated by the plus sign at coordinates (x, y) = (0, 0), and its virial 
radius is marked by the black circle. Also marke d are the v irial radii of the 5 most massive nearby dwarfs and satellites. The metallicity 
is given in units of solar, Zq = 0.0142 (jAsplund et al.H2009l ). 



Here Sij is the trace-free velocity shear tensor and h 
is the measurement scale (taken to be equal to ^sph)- 
This choice for Sij results in no diffu sion for compres - 
sive or purely rotating flows. As in iShen et al.l ((2010), 
we choose a coefficient value of C = 0.05 expected from 
turbulence theory. This diffusion is applied to both 
thermal energy and heavy elements. Metal mixing was 
not implemented in the Eris and Er i sMC simulations of 
iGuedes et all ([201 If ) and IShen etHI ([20121 ). The impact 
of metal diffusion is to spread heavy elements from out- 
flowing wind material to the surroun dings so that m ore 
gas gets contaminated at low levels ([Shen et al.ll2010f ). 
As in previous simulations of this series, the feed- 



back scheme follows the recipes of IStinson etHI ([2006). 
Each Type II SN deposits metals and a net energy 
of 0.8 x 10 51 ergs into a 'blastwave radius', and the 
heated gas has its cooling shut off (to model the ef- 
fect of feedback at unresolved scales) until the end of 
the momentum-conserving phase of the blastwave, which 
is set by the local gas density and t emperature and by 
the t otal amount of energy injected ([McKee fc Ostrikerl 
119771 ). No kinetic energy is explicitly assigned to parti- 
cles within the radius of the blastwave. The energy in- 
jected by many SNe adds up to create larger hot bubbles 
and longer shutoff times. The main difference of this 
feedback model compared to other "sub-grid" schemes 



4 



Circumgalactic medium, galactic outflows, and cold streams 




100 1000 
A [A] 

Fig. 2. The far-UV flux impinging on Eris2's CGM. The 
solid line shows the z = 2.8 UV diffuse background from 
Haardt & Madau (2012), while the colored lines show the local 
stellar radiation at different distances from the center, assuming an 
escape fraction of / e sc = 3% for the ionizing photons that leak from 
the star-forming regions into the CGM. The synthetic galaxy spec- 
trum was produced running Starburst99 (jLeitherer et al. 1 1999) 
assuming a (c onstant) star formation rate of 20 MQyr -1 and a 
IKroiipal OOOlt ) IMF. 

(e.g. Springel fc Hernquistll2003[ ) is that it keeps galactic 
outflows hydrodynamically coupled to the energy injec- 
tion by SNe. In combination with a high gas density 
threshold for star formation (which enables energy de- 
position by SNe within small volumes), this scheme has 
bee n found to be key in p roducing realistic dwarf galax- 
i es ([G overnat o et al. 11201 0| ) and late-type massive spirals 
(jGuedes et al.ll2011[ ). We stress that the parameters of 
the simulation have not been finely tuned to provide a fit 
to any of the observations discussed in this Paper. The 
modifications implemented in Eris2 compared to previ- 
ously published runs of this series have been simply dic- 
tated by improving and updating the physics of galaxy 
formation and the CGM. 

2.1. Local radiation field 

Aside from the isotropic cosmological UV background, 
Eris2's CGM is photoionized by a local non-uniform 
stellar radiation field . We have run Starburst99 
(jLeitherer et al. Ill999[ ) and produced a synthetic galaxy 
stellar spectrum assuming star formation proce eds con- 
tinuo usly at a constant rate of 20 M yr _1 and a lKroupal 
(|2001f ) IMF. The galaxy spectrum, multiplied by the 
distance dilution factor (47rd 2 ) -1 (for d > 0.9 kpc, 
where d is the distance from the center of Eris2) and 
by the frequency-independent (absorption cross-section 
weighted) fraction of ionizing photons that leaks from 
the star-forming regions into the CGM and IGM, (/ eS c)> 
has been added to the UV background to create a non- 
uniform radiation field. This approximation places all 
the stellar sources within the inner kpc, and breaks down 
for distances comparable to the scale length of the stellar 
disk. We have assumed an escape fration of / esc = 3% 
(close to the recent upper lim it, 5%, for z « 3 LBGs de- 
rived by iBoutsia et al. I l2011f ) and run the photoioniza- 
tion code Cloudy to calculate again, in post-processing, 
ionization fractions for an optically thin slab of gas at 
the density, temperature, metallicity, and impinging ra- 



diation field of the simulated SPH particles, under the 
assumption that the temperature of the gas is not mod- 
ified by the local radiation field. 

3. ANALYSIS AND PROPERTIES OF THE CGM 

At z = 2.8 (the lowest redshift reached by our simula- 
tion), Eris2 has a virial mass of M v i r = 2.6 x 10 11 M , 
a virial radius of R v [ r = 50 kpc, a stellar mass of 
M* = 1.5x10^ (a factor 1.4 smaller than in ErisMC, 
see iShen et all l2012f ) and is forming stars at a rate of 
SFR = 20 M yr _1 . The mass loading factor at R v [ r 
(characterizing the amount of material involved in the 
galactic outflow), n = M W /SFR where M w is the rate 
at which mass is ejected, ranges from 0.2 to 0.9 in the 
redshift interval 2.8-3.2. Observations of galactic out- 
flows powered by starbursts suggest a wide range of 
mass loading factors, n = 0.01 — 10, with no obvious 
correlation with the star formation rates of their hosts 
(jVeilleux et al. I [2005). While in low-resolution cosmo- 
logical simulations, the mass loading; factor is one of the 
input parameters (e.g., lOppenheimer fc "D ave 2008]), in 
our SN-driven blastwave feedback scheme there is no spe- 
cific parameter for mass loading. 

Eris2's baryon fraction is 0.13, and the mean gas and 
stellar metallicities within R v { r are (Z g ) = 0.7 Z and 
(Z*) = Z , respectively. The galaxy's H n regions (de- 
fined by all gas at T = 10 4 K within 10 kpc from the 
center) are characterized by 12 + log(0/H) = 8.5, in 
agreement with t he value, 12 + log (O/H) = 8.42 ± 0.06, 
measured by lErb et aT~l (j2006al ) in star- forming galaxies 
of the same stellar mass at (z) = 2.2. About half of all 
the gas-phase heavy elements reside in Eris2's CGM (de- 
fined as all gas at R > 0.2R v { r = 10 kpc), and about half 
of all the CGM metals are locked in a warm- hot compo- 
nent at T > 10 5 K. The total mass of processed material 
with Z > 10~ 3 Z within 150 kpc is 4.7 x 10 10 M , while 
the mass of more pristine, Z < 10 -3 Z , gas within the 
same region is 1.04 x 10 10 M . It is this region that will 
be compared in § HI to observations of the CGM around 
LBGs. Figure [T] shows the projected gas-phase metal- 
licity in a cube 500 (physical) kpc on a side at different 
epochs. While metal-enriched material is seen as far as 
250 kpc from the center, within the same region we iden- 
tify 11 dwarf satellite systems with masses above 10 9 M 
that a re also forming stars and polluting the C GM of 
Eris2 (|Shen et all 120121 : iPorciani fc Madaul l2005f ) . Fig- 
ure [2] shows the far-UV flux impinging on Eris2's CGM. 
With the adopted escape fraction, the galaxy flux at 1 
Ryd dominates over the diffuse UV background at all dis- 
tances < 45 kpc. The most noticeable impact of the stel- 
lar radiation field will be to increase (compared to a case 
with / esc = 0) the interstellar absorption line strengths 
of Si iv and C iv and to decrease the equivalent width of 
hydrogen Lya at small impact parameters. 

The processed simulation provides positions, veloci- 
ties, abundances, and temperatures of the baryonic fluid 
elements represented by SPH particles. To study the 
CGM of Eris2, we draw through the simulation volume 
500 x 500 regularly spaced sightlines within a projected 
sitance of 250 physical kpc from the center, for three 
orthogonal projections of the galaxy. Each particle is 
spread over its 3-dimensional SPH smoothing kernel, and 
a line integral is taken through the smoothed distribu- 
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Fig. 3a. — A map of the projected column density in a cube of 500 (proper) kpc on the side from the Eris2 simulation at z = 2.8. The 6 
panels show H I, C II, C IV, Si II, Si iv, and O VI. Intervals of column density in the range 10 11 - 10 22 cm -2 for H I and 10 11 - 10 16 cm -2 
for all metal ions are marked in the panels with different colors. 
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Fig. 3b. — The multi-phase nature of Eris2's CGM. Same as Fig. [3a] but for cool (T < 10 5 K, top panel) and warm-hot (T > 10 5 K, 
bottom panel) gas only, as traced by the high-ionization species C IV, Si IV, and O VI. 
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Fig. 4b. — Same as Fig. [4a] but for outflowing gas only. 
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Fig. 5. — Distribution of all enriched gas in the temperature-density plane at z = 2.8 within the simulation volume. The top left panel 
shows the mass- weighted total metallicity (in units of solar), while in all the other panels the color coding indicates the ion fractions for 
the species C II, C IV, Si II, Si IV, and O VI. The mean baryon density of the universe at this redshift is 1.4 x 10 -5 cm -3 . Hot enriched 
gas vented out in the halo by the cumulative effect of SN explosions can be seen cooling (first adiabatically then radiatively) and raining 
back onto the disk in a "galactic fountain" . 

tion to determine the line absorption optical depth as 
a function of frequency. For definiteness, the gas-phase 
column density, iV, of an absorbing ion along a sight line 
of length L, is calculated as 



tion can be written as 

r M = ^( m j z j/ m ) w 2B(r j i,h j )a j (u), (4) 



7V = E / (rn 3 Z 3 /m)W(r jh h 3 )dl 
= ^(rnjZ j /m)W 2 B(rji,h j ), 



(3) 



where rrij is the gas particle mass, Zj its mass fraction 
in the relevant ion, and m is the atomic mass of the ion. 
This amounts to a line integral through the 3D smooth- 
ing kernel W(rji,hj) of each particle whose smoothing 
volume Vj = (4/3)7r/i| is pierced by the sightline at 
impact parameter r 3 \. In the equation above we indi- 
cate with W2D the 2D smoothing kernel (with units of 
length -2 ). 

Be cause of the assumed spline smoothing ([Monaghanl 
1992), only SPH particles at impact parameter r 3 \ < 2h 3 
contribute to the line integral in equation (|3]) . Each such 
particle is characterized by a velocity v 3 along the sight- 
line, and by a temperature Tj, so that the total optical 
depth at the observed frequency v in a given line transi- 



and the integrated absorption equivalent width of the 
line is 



-t(v 



Via 



(5) 



Here, the absorption cross-section is expressed in terms 
of the Voigt function as follows: 



Gj(y) 



/+00 
-00 



dy 



exp(-y 2 



(6) 



where s is the frequency integrated absorption cross- 
section (proportional to the oscillator strength / of the 
transition), Avj = uo(bj/c) is the Doppler width, a 3 = 
r/(47rA^j), x 3 = (v — vj) / (Avj) : Vj = v§(\ — Vj/c) is the 
line center frequency corrected for the gas velocity along 
the line of sight, bj = ^/2kTj/m is the Doppler parame- 
ter, and all other symbols have their usual meaning. We 
assume that the galaxy is at zero velocity, and include 
the contribution of bulk motion and Hubble expansion 
to the particle velocity Vj. The flux decrement in equa- 
tion (j5j) is integrated over ±500 kms -1 centered on the 
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Fig. 6. — The covering factor (denned as the percentage of lines 
of sight with a given column density of absorbing gas, as seen by 
the observer) of O VI (top) and C IV (bottom) absorbing material 
through Eris2's CGM as a function of column density N. Values 
are plotted in bins of A log N = 0.1. The color coding marks cov- 
ering factor in different ranges of impact parameter, 0-10 kpc (dark 
blue), 10-50 kpc (cyan), 50-100 kpc (green), 100-200 kpc (yellow). 
Different symbols are used to indicate the median temperature of 
a given column of absorbing gas over all sightlines: T < 10 4 - 5 K 
(squares), 10 4 - 5 < T < 10 5 K (circles), and T > 10 5 K (triangles). 

galaxy's position, and each line of sight passes through a 
simulation volume 1 Mpc (physical) across. 

Figures l3al and l3bl show the total projected column den- 
sities in H i, C ii, C iv, Si n, Si iv, and O vi for Eris2's 
z = 2.8 CGM, and the columns of high ionization species 
that are cool and warm- hot, respectively. Figures Hal and 
l4bl show the same columns for inflowing and outflowing 
CGM material only (gas particles are separated in inflow- 
ing and outflowing according to the sign of their radial 
velocities relative to the center of the main host), while 
the distribution of all enriched gas in the temperature- 
density plane is depicted in Figure El 

Metal-enriched outflows are bipolar and perpendicu- 
lar to the plane of the disk. Heavy elements are clearly 
spread over a large range of phases, from cold star- 
forming material at T < 10 4 K and n > usf — 20 atoms 
cm -3 (corresponding to S = p/p me an ^ 10 615 at z = 2.8) 
to hot T > 10 6 K low density 5 ~ 1 intergalactic gas 
that cannot cool radiatively over a Hubble time. O vi is 
diffuse, has a large covering factor, and extends beyond 
4i? v ir = 200 kpc from the center of Eris2. C iv absorp- 
tion is much clumpier and less extended than O vi, while 



low-ionization species like C n and Si n trace narrow in- 
flowing streams as well as dense outflowing clumps. The 
covering fraction of low-ionization species declines faster 
than that of O vi or C iv. 

Inflows and outflows coexist in Eris2, with about one 
third of all the gas within R v { r found to be outflowing. 
At i?vir, the mean metallicity of inflowing gas is 0.05 Z , 
while the mean metallicity of outflowing gas is 0.56 Zq 
(the average is taken over all gas particles within a thin 
shell of radius R v { r and thickness 0.02 R v [ r ). Inflowing 
enriched gas traces large-scale cold filaments that pene- 
trate deep into the virial radius (see Fig. 2a] and § 14. 2|h 
as well as material raining back onto the disk as part of 
a "galactic fountain". About 77%, 32%, 44%, 66%, 50%, 
and 66% of the total H i, O vi, C iv, C n, Si iv, and 
Si ii mass within 2R Y1T1 respectively, is inflowing, High 
H I column density absorption mainly traces the narrow 
cold inflowing streams. Outflowing material has smaller 
H i columns and contributes to the covering factor of gas 
with A/hi > 10 15,5 cm -2 . 

The distribution of high-ionization species in the 
temperature-density plane is clearly bimodal. At high 
densities and temperatures, ions are in collisional ioniza- 
tion equilibrium (CIE) and their fractional abundances 
are only a function of temperature. At lower densities, 
when photoionization is dominant, the equilibrium abun- 
dances become a function of density as well as temper- 
ature, and metals are typically in a higher ionization 
state at a given temperature than in CIE. C n and Si n 
ions trace primarily cool, 10 4 < T < 10 5 K photoion- 
ized CGM gas at overdensities log S ~2 (as well as dense, 
cold disk neutral material), which extends as far as 2R v [ r 
but with a covering factor that decreases sharply beyond 
R v i r (Fig. [3aj). Si iv and C iv are abundant in low den- 
sity log S > 1 photoionized gas. C iv ions maintain a 
significant covering factor beyond 2R^ 1T and can also be 
detected in warm-hot and denser log S ~ 2 CGM gas 
within i?vir- 

Within 100 kpc, O vi is largely collisionally ionized 
and traces CGM gas with log (5 > 1 and T > 10 5 K 
(in collisional ionization equilibrium, O vi reaches its 
peak abundance fraction, ~ 0.2, at 3 x 10 5 K). A cool, 
T < 5 x 10 4 K clumpy O vi component with smaller 
covering factor becomes dominant at large galactocen- 
tric distances. This is relatively high column, moder- 
ate over density log S ~ 0.5 — 1.5 gas that was enriched 
and expelled by nearby dwarfs and is photoionized by 
the cosmic UV background. The mass ratios of the 
hot versus cold O vi populations are 17:1, 1.6:1, and 
0.7:1 at distances d < i? v ir, ^vir < d < 2R^ 1T1 and 
2ii\dr < d < 3i?vin respectively. 

The covering factor of O vi and C iv absorbing ma- 
terial through Eris2's CGM is plotted in Figure [6] as 
a function of column density. The color coding marks 
the covering factor in different ranges of galactocentric 
impact parameter, while different symbols are used to 
indicate the median temperature of a given column of 
absorbing gas over all sightlines. The figure shows how 
C iv absorbers with large covering factors are typically 
cool (photoionized), while the opposite is true for O vi. 
A hot C iv component is also present at small impact 
parameters and low columns. 

The formalism developed in this Section can be finally 
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Fig. 7a. — A simulated absorption spectrum through Eris2's CGM at z — 2.8, computed at a resolution of 2 kms~ 1 and plotted in 
velocity space. A random sightline at impact parameters b = 0.46i? v i r = 23 kpc from the galaxy center is shown. Gas with positive radial 
velocity with respect to systemic is redshifted. The total column densities of the absorbing H I, C II, C IV, Si II, Si IV, and O VI ions are, 
respectively 2.6 x 10 16 ,1.1 x 10 14 ,4.2 x 10 14 ,8.7x 10 12 ,5.3 x 10 13 , and 1.3 x 10 14 cm" 2 . 
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Fig. 7b. — Same as Fig. I7al but for a sightline at impact parameter b = 34 kpc. The total column densities of the absorbing H I, C II, 
C IV, Si II, Si IV, and O VI ions are, respectively 1.3 x 10 17 ,4.2 x 10 14 ,4.6 x 10 i4 ,3.7 x 10 13 ,1.3 x 10 14 , and 1.8 x 10 14 cm- 2 . 
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Fig. 8. — Total projected metal mass density, along each 

sightline through Eris2's CGM, as a function of galactocentric im- 
pact parameter, 6, in units of the virial radius. The metal mass 
density is plotted at two different redshifts, z = 2.8 (i? v ir — 50 
kpc) and z = 5 (i? v ir — 19 kpc). The density of black points (one 
for each sightline) is a measurement of the covering factor of heavy 
elements in the simulation. The solid yellow line shows the median 
value of at a given b/R v - 1Y . At the two redshifts, the metal 
bubbles around Eris2 appear scaled versions of each other. 



used to compute the optical depth of the hydrogen Lya 
A1216, C ii A1334, C iv AA1548,1550, Si n A1260, Si iv 
A1393 and O vi A1032 transitions. Figures EH and l7bl 
show artificial absorption spectra along two lines of sight 
through the CGM at impact parameters b = 23 kpc 
and b = 34 kpc from Eris2's center, plotted on a galac- 
tocentric velocity scale. Column densities range from 
8.7 x 10 12 cm -2 for Si n to 1.3 x 10 17 cm -2 for hydrogen 
Lya, and individual components span a radial velocity 
range from about —200 to +200 kms -1 . Velocity profiles 
are complex, with ouflowing and inflowing gas detected 
at positive (redshifted) and negative (blueshifted) pro- 
jected velocities from the near and the far side of the 
galaxy. The simultaneous presence, in velocity space, 
of features from singly ionized species as C n and Si n 
along with higher-ionization species reflects the multi- 
phase structure of the simulated CGM. The sightline de- 
picted in Figure [7b] goes through cool, metal-enriched 
gas infalling on the near side of the galaxy. The inflow- 
ing component is distributed in physical space between 
R v i r and 2R^ 1T1 has overdensities log 5 > 2, metallicities 
in excess of 0.03 Z@, and is part of the metal-enriched 
zone around a dwarf satellite companion accreting along 
a filament. We will discuss the observability of cold acc- 
cretion streams via low-ionization metal lines in the next 
section. The spectrum shows several outflowing com- 
ponents, including solar-met allicity photoionized gas at 
v ~ —150 kms -1 that is being ejected by the main host, 
as well as pristine IGM gas at v ~ kms -1 that is lo- 
cated > 240 kpc from the galaxy center. 

A detailed analysis of interstellar absorption line 
strengths as a function of redshift will be the subject 
of another paper. Here, it is important to keep in mind 
that, when properly scaled, the property of Eris'2 CGM 
are not changing rapidly with cosmic time. Figure [8] 
shows the total projected metal mass density, along 



each sightline, as a function of galactocentric impact pa- 
rameter in units of the virial radius, b/R v { r . The metal 
mass density is plotted at two different redshifts, z = 2.8 
(R v [ r = 50 kpc) and z = 5 (i? v ir = 19 kpc). The den- 
sity of black points (one for each sightline) is a measure- 
ment of the covering factor of metal-enriched material in 
the simulation. At the two redshifts, the metal bubbles 
around Eris2 appear scaled versions of each other. 

4. COMPARISON WITH THE OBSERVATIONS 

4.1. Circumgalactic metals at z=2-3 

The close pair s ample of galaxies at (z) = 2.2 of 
iSteidel et al.l (|201Q[ ). drawn from a spectroscopic survey 
of "BX" UV-color selected objects, provides a robust 
map of cool circumgalactic gas at galactocentric impact 
parameters b = 3 — 125 kpc. Objects in this sample have 
been well-characterized in terms of the ir spatial cluster- 
ing and stellar population parameters ([Adelberger et all 
l2QQ5al lErb et al. ([2006bD . A good match to both the 
observed clustering strength and the space density of 
the "BX" galaxies to dark matter halos in the Mille- 
nium simulation is obtain ed for an ave rage halo mass of 
(M vir ) = 9 x 10 11 M (jConrov et al. Il200l . The virial 
radius for a halo at the average mass that formed sev- 
eral 10 8 yr prior to the epoch of observation is R v [ r = 80 
kpc (Steidel et al. 2010), close to the scale at which low- 
ionization absorption lines become weaker. This sug- 
gests that beyond R v { r absorbers may become too dif- 
fuse and highly ionize d to produce signific ant columns of 
low-ionization metals (St eidel et al.|[20To[ ). With a virial 
mass at z = 2.8 of M Y1Y = 2.6 x 10 11 M and a ra- 
dius i^yjr = 50 kp c, Eris2 is a smaller counterpart to 
ISteidel et all (|2010f ) close pair sample galaxies |j For a 
proper comparison between simulation data and obser- 
vations, we have therefore scaled the impact parameter 
with the virial radius of the mean halo in the pair sam- 
ple (80 kpc) and with the Eris2's virial radius in the 
simulation, and plotted the observed and predicted line 
strengths against the correspondingly normalized impact 
parameters. Along every simulated sightline, the rest- 
frame equivalent widths Wo of C n 1334, C iv 1549 (in- 
cluding both components of the doublet), Si n 1260, and 
Si iv 1393 are shown in Figure [9] versus b/R v [ r: together 
with the values observed in the composite spectra of 
foreground galaxies (placed at (b) = 1 kp c) and back- 
groun d galaxies (at (b) = 31, 63, 103 kpc) by St eidel et al.l 
(|2010l ) . The density of black points is a measurement of 
the covering factor of absorbing material at a given im- 
pact parameter in the simulation. We plot Wo along ev- 
ery artificial line of sight rather than an average over all 
sightlines in order to underline the dispersion in the mea- 
sured line strengths as well as the contribution of satellite 
galaxies (the narrow "spikes" in equivalent width seen at 
b/Rvir > 1)- For each species, the strength of absorp- 
tion declines rapidly around b =1-2 i? vir . The agreement 
between the simulation and the observations is rather 
good: except perhaps for a slight deficit in the covering 
factor of C n and Si n around ~ 0.8R v { r and an excess of 
C iv at large impact parameters , Eris2's simulated CGM 
broadly reproduces the line strengths-impact parameter 

6 This is also true in the luminous component, as the mean stellar 
mass of the spectr oscopic sample, (M*) = (3.6 ± 0.4) x 10 10 Mq 
(|Erb et al. 112006131 ). is 2.5 times larger than in Eris2. 
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Fig. 9. Interstellar absorption line strengths of the 1216 Lya, 1334 C II, 1549 C IV, 1260 Si II, 1393 Si IV, and 1032 O VI transitions 
as a function of galactocentric impact parameter b in physi cal kpc. Valu e s for C IV include both components of the AA1548,1550 doublet. 
The cyan points at 6=1, 31, 63, and 103 kpc are data from Steidel et al. (2010), with downward arrows on points indicating upper limits. 
The yellow points in the top left panel are data from Rakic et al. (2012). The impact parameter has been normalized to R V1Y = 80 kpc 
(see text for details). The black dots show the rest-frame equivalent widths measured along artificial sightlines through the CGM of our 
Eris2 simulation at z = 2.8, and the yellow lines marks the median and the 25th and 75th percentiles. 



trends observed by Steid el et al.l (|2Q1Q[ ). In particular, 
the covering factor of absorbing material declines less 
rapidly with impact parameter for Lya and C iv (and 
O vi) compared to C n, Si iv, and Si n, with Lya re- 
maining strong (W^ya > 300 mA) to > 5R v [ r = 250 kpc. 
At small impact parameters, where lines are mostly sat- 
urated, the equivalent width along a sightline is mod- 
ulated by the velocity structure of the absorbing mate- 
rial. The kinematics of Eris2's inner CGM appears there- 
fore to be consistent with the observat ions. A compari- 
son be tween Figure [9] and Figure 15 of iFumagalli et al~l 
(|2011l ) shows that the predicted interstellar Si n equiv- 
alent widths are about an order of magnitude larger in 
Eris2 than in IFumagalli et alTI (|2011f ). and therefore that 
strong galactic outflows are an essential ingredient for a 
quantitative modeling of the CGM. 

4.2. Circumg atactic H i and cold flows 

Figure [9] also shows the equivalent width of hydro- 
gen Lya absorption as a function of impact parameter 
for the galaxy-galaxy pair sample of iSteidel et al.l ([2010) 



as we ll as the QSO-galaxy pair sample of iRakic et al.l 
(2012]). A gain, the simulation agrees well with the ob- 
servations, as Lya remains strong (W^ya > 300 mA) 
to 5i?vir = 250 kpc. Infall of cool gas via nar- 
row streams, directly onto the central regions of galax- 
ies, is the mode of gas accretion that is predicted 
to feed st ar formation in h i gh redshift galaxies like 
Eris2 (e.g. iKeres et al.l |2QQ5t iDekeTfc BirnboIm| 120061: 
IQcvirk et an 120081 : lAgertz. Tevssier. fc Moorel 120091 : 
iDekel et al. 1120091 ). Observational signatures of cold ac- 
cretion streams at high redshift are key for testing our 
understanding of the flows of matter, energy, and metals 
into and out of galaxies. Absorption by neutral hydrogen 
(and, as we shall argue below, by low-ionization metal 
lines) is a promising way to observe the cold streams. 
Recent observations of 2.0 < z < 2.8 star- forming galax- 
ies near a QSO sightline (the Keck Baryonic Structure 
Survey) show that the covering fraction of absorbers of 
various A^hi columns around galaxies ranges from 90 ±9% 
(< i? vir ) to 68 ± 9% (< 2R vir ) for logiV H i > 15.5, and 
from 30 ± 14% (< R vir ) to 28 ± 9% (< 2R vir ) for Lyman 
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Fig. 10. — The covering factor at different galactocentric impact 
parameters of absorbers of various species and column densities (in 
cm" 2 ) in Eris2's CGM. 

Limit Systems (LLSs) with logA^i > 17.2 (jRudie et al.l 
l2012h . 

The simulation (see Fig. [TQ|) yields covering factors of 
71% and 27% (< R vir ) and of 41% and 12% (< 2R vir ) 
for log iVm > 15.5 and > 17.2, again in reasonable agree- 
ment with the data, albeit perhaps on the low side for 
the la rgest impact parameters. According to lRudie et al.l 
(2012), the incidence of absorbers with logA^i > 14.5 
is found to be higher near galaxies than in the gen- 
eral IGM. The covering factor of log 7Vhi > 14 systems 
within 2.5 R v [ r is observed to be 92 ± 5% vs. 99% in 
Eris2. Eris2's covering factor of LLSs appears compara- 
ble or higher than those predicted at high redshift by 
IFaucher-Giguere fc Keresl (|2Qlll ) and iFumagalli et al I 
(|2Qllf ) for galaxies in the same mass range. Both these 
groups ran high-resolution cosmological hydrodynami- 
cal si mulations that do not g enerate strong galactic ou- 
flows ([Fu magalli et a l7]l2Qlll u sing the AMR technique, 
IFaucher-Giguere fc Keresl I2Q11I using SPH as in Eris2), 
combined with radiative transfer. Because our hydrody- 
namical simulations lack proper ionizing radiative trans- 
fer, we cannot correctly address the covering factor of 
the self-shielded, thicker Damped Lya (DLA) absorbers 
with log 7V H i > 20.3. 

It is interesting at this stage to look at the properties 
of cold streams. About 90% of all optically thick absorb- 
ing gas within 1-2 virial radii is found to be inflowing 
wit h rad ial velocities < 150 — 200 kins" 1 . Figures Tllal 
and lllbl shows two projections of all inflowing material 
with logA^Hi > 17.2, as traced by H i and C n, respec- 
tively. The streams can be clearly seen penetrating deep 
inside the virial radius and delivering cold gas to the cen- 
tral galaxy. Enriched to metallicities above 0.01 solar by 
previous episo des of star format i on in the main host an d 
nearby dwarfs (|Shen et al.l 120121 : IFumagalli et al. Il2011f ). 
the streams give origin to strong (Nqu > 10 13 cm -2 ) 
C ii absorption with a covering factor of 22% within 
R v i r and 10% within 2i? v i r , which should make the pres- 
ence of cold flows detectable with metal absorption lines. 



We find no substantial suppression of the cold accre- 
tion mode caused by galactic-scale mass outflows. The 
mass inflow rate of cold (T < 10 5 K) gas into Eris2 is 
^coid = 18 M^yr~ * (measured at i? v ir), in agreement 
with the results of iKeres et al. I (j2009f ) for galaxies in 
the same redshift and mass range as Eris2, and compa- 
rable to Eris2's star formation rate. About 35% of the 
cold gas mass is brought in by dwarf satellite galaxies. 
The cold mode dominates over the hot mode, which is 
characterized by an accretion rate of Mhot = 5 M yr _1 . 
The distribution of metallicity, column density weighted, 
in the entire sample of inflowing LLSs within 2R v [ r is 
centered around a median value of 0.07 Z^, significantly 
higher than found in t he sim ulations of IFumagalli et al71 
(|2011f ): iKimm et all ()2Qllh . Nearly all inflowing LLSs 
have metallicities above 0.01 solar (Figure [T2]). 

4.3. The vi halos of star-forming galaxies 

As already discussed in the previous section, as a 
consequence of the SN blastwave- driven galactic out- 
flows generated in our simulations, the central galaxy 
is surrounded by a large halo of collisionally ion- 
ized O vi. Ubiquitous, large (150 kpc) halos of O vi 
around low-redshift star-forming galaxies have been re- 
cently detected in a survey with th e Cosmic Origins 
Spect rograph onboard of the HST (ITumlinson et al.l 
120111 ). and it is interesting t o compare their proper - 
ties with our simulation. The ITumlinson et al. I (|2Qllh 
sample includes galaxies with stellar masses 9.5 < 
log(M*/M ) < 11.2 and specific star formation rates 
-11 < log (sSFR/yr _1 ) < —9. Eris2's stellar mass and 
sSFR at z = 2.8 are in the middle and on the high 
side of such ranges, respectively. We have computed the 
total O vi column along simulated sightlines at differ- 
ent impact parameters a nd plotted the results in Fig- 
ure [T3l together with t he ITumlinson et al.l (j2011f ) and 
iProchaska et al. I (|2011l ) low-redshift observations. Our 
high-redshift simulated galaxy again matches the trends 
observed in similar systems at low redshfts, with typical 
column densities of A/ovi ^ 10 14 cm -2 and a near unity 
covering factor are maintained all the way out to 100-150 
kpc. 

5. SUMMARY 

Recent observations of the CGM of star- forming galax- 
ies at high and low redshift have called for more strin- 
gent comparisons with hydrodynamical simulations of 
galaxy formation. The high-resolution "zoom-in" sim- 
ulation presented here, Eris2, adopts a blastwave scheme 
for supernova feedback which, in combination with a high 
gas density threshold for star formation, has been shown 
to generate large-scale galactic outflows and be key in 
producing realist ic dwarf galaxies and l ate-type mas- 
sive spirals (e.g. iGovernato et al.l 120101 : iGuedes et al.l 
20111). in turning dark m atter cusps into cores (e.g. 
Pontzen fc Governatoll2012[). and in enr iching the CGM 
of high redshift LBGs (|Shen et al.ll2012T ). Simulations of 
star-forming galaxies run to z = have recently stressed 
the importance of strong stellar feedback in reproduc- 
ing the O vi hot gas eous halos observed in local galaxies 
(jStinson et al.ll201l| ). Eris2 includes metal-dependent ra- 
diative cooling, a model for the diffusion of metals and 
thermal energy, and a local UV radiation field added in 
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Fig. 11a. — Cold accretion streams in Eris2: a map of hydrogen column density showing inflowing, optically thick material with 
logA^Hi > 17.2 for 2 different projections. Eris2's virial radius is marked by the white circle. 





Fig. 11b, 

post-processing, and is ideally suited for a detailed study 
of the kinematics, thermal and ionization state, and spa- 
tial distribution of metal-enriched gas in the CGM of 
massive galaxies at redshift ~ 3. The results presented in 
this paper do not depend sensitively on metal diffusion: 
a twin simulation without diffusion produces a slightly 
clumpier CGM, but the impact on the covering factor of 
metal ions is small. 

We have generated synthetic spectra by drawing sight- 
lines through the simulated CGM at different galactocen- 
tric impact parameters, and compared the theoretical in- 



terstellar absorption line strengths with the observations. 
We emphasize again that the parameters of the simula- 
tion have not been tuned to provide a fit to observations 
of the CGM. In our examination we have found that: 

• Eris2's CGM at z = 2.8 contains multiple phases 
having a wide range of physical conditions, with 
half of all gas-phase metals locked in a warm-hot 
component at T > 10 5 K. Outflows and inflows co- 
exist, and about one third of all the gas within 
R v [ r is outflowing. Inflows bring in material along 
optically thick "cold" streams. The streams are en- 
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Fig. 12. — Cold, optically thick streams in Eris2 are metal- 
enriched. The cumulative covering factor of logA^Hi > 17.2 in- 
flowing material above a metallicity Z is plotted within R V1T and 
2R vir . Streams within R v [ r (2R V - 1Y ) are enriched above 0.03 (0.01) 
solar by previous episodes of star formation in the main host and 
in dwarf satellite galaxies. 

riched to metallicities above 0.01 solar by previous 
episodes of star formation in the main host and in 
nearby dwarfs. 

• The C ii and Si n ions trace primarily cool, 10 4 < 
T < 10 5 K photoionized CGM gas at overdensities 
log S ~2 extending as far as 2i? v i r , but with a cover- 
ing factor that decreases sharply beyond i? v ir- The 
distributions of high-ionization species are clearly 
bimodal, with ions in CIE at high densities and 
temperatures, and with photoionization becoming 
dominant at low densities. C iv ions maintain a sig- 
nificant covering factor beyond 2R Ylv and can also 
be detected in warm- hot and denser log 6 ~ 2 CGM 
gas within i? v ir- The O vi diffuse halo extends be- 
yond 4i? v i r = 200 kpc from the center. Within 
100 kpc, O vi is largely collisionally ionized and 
traces CGM gas with log 6 > 1 and T > 10 5 K. 
A cool, clumpy O vi component, photoionized by 
the cosmic UV background, becomes increasingly 
important at larger distances. 

• Synthetic spectra generated by drawing sightlines 
through the simulated CGM produce interstellar 
absorption line strengths of Lya, C n, C iv, Si n, 
and Si iv as a function of galactocentric impact 
parameter (scaled to the virial radius) that are 
in good agreement with those obser ved at high- 
redshift by Steidel et al. (2010) (and iRakic et al.1 
2012 in the case of Lya). The covering factor of 
absorbing material is found to decline less rapidly 
with impact parameter for Lya and C iv (and 
O vi) compared to C n, Si iv, and Si n. Lya re- 
mains strong to > 5R v { r = 250 kpc. At small im- 
pact parameters, where lines are mostly saturated, 
the equivalent width along a sightline is modulated 
by the velocity structure of the absorbing mate- 
rial. The kinematics of Eris2's inner CGM appears 
therefore to be consistent with the observations. 
Note that, contrary to the spherically outflowing 



wind model of iSteidel etHI (j2010f ). in our simula- 
tions inflowing material contributes significantly to 
the absorption in the Wo — b plane. 

• At small impact parameters, the local galactic 
UV flux increases the abundance of high-ionization 
species like C iv and Si iv, and decreases that of 
C ii, Si ii, and Hi. O vi is largely collisionally 
ionized within 100 kpc, so the effect of the local 
UV on O vi is small. More quantitatively, within 
i? v i r , galactic radiation increases the total amount 
of O vi by 1%, C iv by 128%, and Si iv by 146%, 
while it decreases C n by 53%, Si n by 51%, and H i 
by 87%. Nevertheless, the effect on the equivalent 
widths of Lya and metal transitions is small. This 
is because absorption lines are saturated in the lo- 
cal CGM and their equivalent widths are mainly 
determined by the kinematics and covering factor 
of the gas rather than by its column density. 

• The fraction of sightlines within one virial radius 
that intercept optically thick, TVhi > 10 17,2 cm -2 
gas is 27% , in good ag r eemen t with recent obser- 
vations bv lRudie et"aLl (|2012h who find 30 ± 14%. 
Within 2R v [ r , our co vering factor i s 12% , below the 
28 ± 9% observed bv lRudie et al.l (|2012[ ). The cor- 
responding covering factors for logA^Hi > 15.5 are 
(71%, 41%) versus (90 ± 9%, 68 ± 9%) found in the 
observations. 

• Optically thick material is found to trace cold 
accretion streams that penetrate deep inside the 
virial radius and deliver cold gas to the central 
galaxy. The streams have metallicities above 0.01 
solar. C ii absorption with iVcn > 10 13 cm -2 arises 
from such streams with a covering factor of 22% 
and 10% within one and two virial radii, respec- 
tively, making the presence of cold flows detectable 
with metal absorption lines. There is no substan- 
tial suppression of the cold accretion mode caused 
by galactic outflows. 

• The O vi halo gives origin to CGM absorption with 
a typical column density of iVovi ^ 10 14 cm -2 
and a near unity covering factor that is main- 
tained all the way out to 150 kpc. This matches 
the trends recently o bserved in star-forming g alax- 
ies at low redshif t bvlTumlinson et al. I (|2011[ ) and 
iProchaska et al. I (|2011f ). 

While our zoom-in simulations of this single system ap- 
pear then to reproduce quantitatively the complex bary- 
onic processes that determine the exchange of matter, 
energy, and metals between galaxies and their surround- 
ings, we acknowledge a number of caveats to the work 
presented here. Radiation transport is not explicitly 
present in our simulations, and while we have checked 
the robustness of our conclusion with a self-shielding ap- 
proximation used in post-processing, simulations that in- 
corporate some form of radiative transfer are in the mak- 
ing. Eris2's interstellar absorption features do not reach 
the maxi mum velocitie s \v m ^ I — 800 kms -1 observed in 
LBGs bv ISteidel et al.l J2010), and zoom-in simulations 
of more massive systems are needed in order to fully test 
our galactic outflow model against observations of the 
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Fig. 13. — The O VI column density vs. impact paramete r for the star-formin g subsample of Tumlinson et al. (2011) (cyan squares ) and 
the L * (red dots) an d sub-L* (red triangles) s ubsamples of Prochaska et al. (201J]) of low-redshift galaxies. The L* data of Prochaska et al. 
(2011) and those of Tumlins on et al.1 {2011) have been scaled assuming a virial radius of 250 kpc. The sub-L* data of Prochaska et al. 
(2011) have been scaled assuming a virial radius of 160 kpc. The black dots show the O VI column density measured along artificial 
sightlines through the CGM of our Eris2 simulation at z = 2.8, and the yellow lines marks the median and the 25th and 75th percentiles. 



high-redshift CGM. The dumpiness of Eris2's CGM may, 
at least partly, be a numerical artifact of the SPH tech- 
nique. A clumpy, two-phase medium could be produced 
by thermal in stabilities in the warm/hot halo around 
the main ho st (|Maller fc Bullo eld 12003 : iKaufmann et al. I 
2006, 120091 ), although whether or not the conditions for 
thermal instabilities to op erate are met in real galaxies is 
still a controversial issue (jBinnev et al. 1 12009). The pro- 
cess has been observed in idealized simulations in which 
gas cools out of an initially uniformly hot corona and the 
instability is seeded by Poisson-noise density fluctuations 
i n the initial particle distributions. iKeres fc Hernquistl 
(2009) have argued that clouds arise naturally in realis- 
tic cosmological hydrodynamical simulations as a result 
of a combination of Ray leigh- Taylor instabilities at the 
interface between cold accretion flows and the hot phase 
and thermal instabilities. Moreover, hot outflows and 
winds from supernovae launched from the disk might in- 
duce strong, large wavelength fluctuations that are not 



easily damped, and may thus act as the seeds to thermal 
instability in the CGM. These seeding mechanisms, as 
well as perturbations by satellites, may all play a role in 
our simulations. A detailed analysis of the dumpiness 
of Eris2' CGM as the possible result of such instabilities 
will be the topic of future work. 
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